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Materials Modeling at NASA ARC 


New materials modeling: interdisciplinary 

• Computational chemistry 

• Computational physics 

• Computational engineering 
Ablative composites 

• Application: atmospheric re-entries 

• Materials: PICA, Avcoat,... 

Ultra high temperature ceramics (UHTC) 

• Applications: leading edges 

• Materials: ZrB 2 and HfB 2 

UHTC modeling illustrate our approach 



PICA: fibers w/ phenolic 



UHTC: HfB 2 w/ 20% SiC 




UHTC for Sharp Leading Edges 


Sharp leading edge for hypersonic aircraft 

- Enhances vehicle performance 

- Improves safety 

Higher temperature requirements 

- Shuttle RCC leading edge: T~1650C 

- Sharp leading edged vehicles: T>2000C 
UHTC advantages for sharp leading edges 

- Reasonable mechanical properties 

- Oxidation resistance 

- Hjgh themial conductivity 

- Effective thermal radiation 

- Thermal shock resistance 



Leading edges of hypersonic vehicle 



SEM image of UHTC microstructure 


Gasch, Johnson and Marschall, J. Am. Ceram. Soc. 91, 1423, (2008) 
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Multiscale Modeling of UHTC 



* Materials modeling will: 

- Discover fundamental mechanisms 

- Structure-property relationships 
— Design new materials 

- Accelerate material development 

* Framework integrates three methods: 

- Ab initio - fundamental chemistry 

- Atomistic- thermal/mechanical 

- Continuum - macro properties 

* This talk focuses on thermal conductivity 

- Atomic structure and bonding 

- Interatomic potential development 

— Lattice thermal conductivity simulations 

- Grain boundary thermal resistance 

- Imaged based FEM of GB networks 




!•••••! 

• • • • • 




JL, Daw, Squire and Bauschlicher, (2012), submitted 
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NASA 




UHTC: ZrB 2 









Alternating layers of 
Zr/Hf (red) and Boron (gray) 


and HfB 


r*Y#Y*Y«r*Y* 

r^T^rir^T^ 


Graphitic Boron layers 
with Zr/Hf over each ring 



NASA 


Electron Localization Function 





Covalent bonding in Boron plane Metallic bonding in Zr planes 


Blue = High 
Red = Low 


Ionic bonding between Interlayers 


JL, Bauschlicher, and Daw, J. Am. Ceram. Soc, (2011) 



Fundamental Properties: ZrB 2 & HfB 2 



Electronic Spectra 



E-Ef(eV) 

Electronic properties 
essentially identical 


Vibrational Spectra 



Frequency (cm-1) 

Vibrational differences due 
to Zr/Hf mass difference 



Tersoff Bond Order Potential 


• Two body terms (A,X, B, jn) energy 

/•- Zl// 

i*j 

f^(d) = A ij e^p(-A iJ d) 
fA ] (d) = -B ij exp(-v ij d ) 


• Bond order (P, ^ 3 , n, m) 

fy=c :i+A"'Q)~ /2ni 

iij = Z fc Jl ( r ik ^YijkSi Wqk ) eX P [^ 3 , (djj ~ d ik )”' ] 

k*i,j 

• Angular function (c, d, h) 

g t (0)= 1 + cf / df - cf /[df + (h t - COS 0 2 ] 



Interatomic Energy 


Daw, JL and Bauschlicher, Comp. Mat. Sci., (2011) 
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Lattice Thermal Conductivity 


• Green-Kubo thermal conductivity tensor 

B 


• Heat current J (x^vj, energy e t , stress-tensor S, 



X e > v > ~H s i v i 

i i _ 

i ^ i<j 


JL, Daw and Bauschlicher, J. App. Phys, (2011) 
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<J(t)J(0)>/<J(0)J(0)> 



Heat Current Correlation Function 
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HfB; 300K 
liZrB ? 1000K 
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15 


20 


• Monoatomic systems (e.g. Si) have monoatomic decay 

• ZrB 2 has longer period than HfB 2 at T=300K 

• ZrB 2 at T=1000K has longer period than T=300K 


JL, Daw and Bauschlicher, J. App. Phys, (2011) 


10 


Power (ar bitary units) 


Correlation Function Power Spectra 



• Correlations oscillates with metal-B optical modes 

• C xx and C yy oscillate with in-plane mode frequency 

• C zz oscillates with out-of-plane mode frequency 
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K(W/(m.K)) 


Lattice Thermal Conductivity: ZrB 



• 8 independent, 10 ns simulations, T=300K 

• 8x8x16 unit cell, 12,255 atoms 

• k x =60 W/(m.K), k zz =40 W/(m.K) 


K(W/(m.K)) 


Lattice Thermal Conductivity: HfB 



• 8 independent, 10 ns simulations 

• 8x8x16 unit cell (12 atoms)= 12,255 atoms 

• k- xx =76 W/(m.K), k z =65 W/(m.K) 


k (W/(m.K)) 


Thermal Conductivity vs Temperature 



•8 independent, 10 ns simulations for each point 

• Data fit to 1/T curves 
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Single Crystal ZrB 2 Data 


• Electron & lattice thermal conductivity 


K — K + K 


lat 


• K| at is 0.3k from polycrystalline data 

• Single crystal data 

• k xx = 140 W/mK, k zz = 100 W/mK 

• 1 sample, 1 measurement 

• Defects uncharacterized 

• k xx = 45 W/mK , k, z = 30 W/mK (lattice) 

• More data needed for ZrB 2 and HfB 2 


Kinoshita, Otani, Kamiyama, Amano, Akasaki, Suda, Matsunami, Japan. J. App. Phys v (2001) 
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UHTC Grain Boundaries 



Z7 symmetric tilt 
(graphene GB structure) 



YJ symmetric twist 


Full ab initio / MD analysis of two tilt and two twist boundaries 



JL, Daw, Squire and Bauschlicher, J. Am. Ceram. Soc, submitted, (2012) 
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Interface Thermal Resistance 


As +As -As 


Jz 

1 s 

Jz 

Jz 

Jz 

N. 1 






GB1 GB2 

Swap atoms to create heat source/sink 


Q = <j k A T 

g k is the Kapitza conductance 


Temperature Profile 



GB at z = 0.25 and 0.75 


Muller-Plathe, J. Chem. Phys. 106 (1997), p 6082 
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Simulation Results 



Interface Conductance 





Ab Initio 

DLB/Pot 



7 (meV/A 2 ) Az (A) 

7 (meV/A 2 ) A* (A) 

ok (GW/(m 2 -K}) 

c-tilt 

153(309) 

112(233) 

1.77 

c- twist 

157(375) 0.20 

111(258) 0.29 

0.58 

a- tilt 

227(1040) 

107(1380} 

0.3S/0.31 

a- twist 

212(1230) 

118(1430) 

0.55/0.53 




TABLE I: Energetics and thermal conductance for ZrB^ grain hcmna arv at r ue lu re s from 
empirical potentials (DLB/Pot 1) and ah mitjo/DFr. Units for 7 are meV/A 3 , A z are A 
and ur are GW/{m 2 - K). 


• Very high thermal conductance (very low resistance) 

• Experimental data indicate much lower values 

• Not surprising given pristine grain boundaries 

• Need: improved processing for improved GB properties 

• Need: modeling more complex boundaries 


Experimental Results: Polycrystalline ZrB 2 



ZlJft 400 6«fJ SCO KKK) I20U NON 1«(W ISOO 

Temperature (K) 

At 300K, K tot = 55 W/mK, K e = 33 W/mK, K lat = 22 W/mK 


Fahrenholtz, Hilmas, Talmy and Zaykoski, J. Am. Ceram. Soc., (2007) 

Zimmermann, Hilmas, Fahrenholtz Dinwiddie, Porter, Wang , J. Am. Ceram. Soc., (2008) 
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Microstructural Model: ZrB 2 



What is effect of grain boundary network 
on thermal conductivity? 


Estimate with Brick Layer Model 




B1 B2 B3 BN 




AT 0 d 


K eff = 48 W/mK 
using 

Ko= 50 W/mK, 

R K = 1 m 2 K/GW, 
d = 6 jam 


Very small 
reduction using 
MD resistances 
and BLM! 


Fahrenholtz, Hilmas, Talmy and Zaykoski, J. Am. Ceram. Soc., (2007) 
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Imaged based FEM 




boundaries 


< 51.6(j,m ► 

• Larger reduction with realistic structures and parameters? 

• Realistic microstructure 

• Finite element mesh from OOF2 

• MD thermal conductivity for grains 

• Experimental interface resistance 


Reid, Lua, Garcia, Coffman, Langer, Int. J. Mat. and Prod. Tech., 2009 
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Development of Steady State 



Inc: 0 

Time: 0.000e+00 


MS^, Software 


0 


301.0 
300.9 

300.8 
300.7 
300.6 
300.5 
300.4 
300.3 
300.2 

300.1 
300.0 

299.9 



ABAQUS-STYLEFI LECREATEDBY00F20N2010-09-1318 : 08:15. 781461FR0M 
Temperature 1 


Uniform thermal gradient (UTG) applied vertically across structure 
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Effective Thermal Conductivity 


• Boundary conditions 

• Uniform temperature gradient (UGT) 

• Uniform heat flux (UHF) 

• Transport direction 

• Vertical 

• Horizontal 

• Evaluate effective properties 


<«}=-v< vr ) 


k eff = effective thermal conductivity 
(q) = volume averaged heat flux 
(Vr) = volume averaged temperature gradient 

• Is microstructure "representative"??? 


k uhf — k — k utg 





301.0 
300.9 
300.8 
300.7 
300.6 
300.5 
300.4 
300.3 
300.2 

300.1 
300.0 








Effective Thermal Conductivity 



BC Type 

Parameter 

Vertical 

Horizontal 


<q> 

- 27.32 

19.05 

UGT 

<dT/ds> 

1.56 

- 1.17 


K eff 

17.48 

16.24 


<q> 

- 28.13 

28.05 

UHF 

<dT/ds> 

1.68 

- 1.76 


K eff 

16.72 

15.93 


Comparison of results: 

• BLM K eff = 15.4 (series model) 

• Rule of mixtures K eff = 44.14 (parallel model) 

• FEM has series and parallel contributions 

• BLM has very good agreement with FEM 
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Conclusions 



• NASA ARC computational materials modeling: 

• Ablative composites 

• Ultra high temperature ceramics 

• Multiscale framework for UHTC: 

• Ab Initio - bonding, electronic & vibrational spectra 

• Atomistic simulation - bulk and interfacial thermal conductivity 

• Continuum - microstructural modeling and effective properties 

• Iteration with experiment needed to "close" loop 

• Modeling unanswered questions: 

• Interatomic potential fidelity 

• Complex grain boundary structural models and properties 

• Experimental unanswered questions: 

• Single crystal thermal conductivity 

• Electronic vs lattice carrier breakdown 

• Grain boundary atomic structures and properties 

• Improved grain boundaries from improved processing 


25 


